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Abstract 

We study effects of cosmic fluids on finite-time future singularities in modified f(R, G)- 
gravity, where R and G are the Ricci scalar and the Gauss-Bonnet invariant, respectively. 
We consider the fluid equation of state in the general form, uj — u)(p), and we suppose the 
existence of a bulk viscosity. We investigate quintessence region (uj > —1) and phantom region 
(oj < —1) and the possibility to change or avoid the singularities in f(R, G)-gravity. Finally, 
we study the inclusion of quantum effects in large curvature regime. 

I Introduction 

Recent observational data imply that the current expansion of the universe is accelerating [J [5]. 
This is the so called Dark Energy issue. There exist several descriptions of this fact. Among them, 
the simplest one is the introduction of small positive Cosmological Constant in the framework of 
General Relativity (GR), the so called A-CDM model. A generalization of this simple modification 
of GR consists in considering modified gravitational theories, in which the action is described by 
the Ricci scalar R plus an arbitrary function f(R) of R (for reviews, see [3]). 

Alternatively, accelerating FRW universe may be described by cosmological quintessence/ 
phantom scalar dark energy or GR plus fluid, satisfying suitable equation of state. It is also 
well-known that any of such dark energy-models maybe represented as an effective fluid with 
corresponding characteristics. 

In this paper, we would like to investigate modified f(R, G) theories of gravity, where G is the 
Gauss-Bonnet invariant. These models are generalization of f(R) models, and they mainly are 
inspired by (super) string theories [HE]. 

It is also well-know that many of such modified gravity models bring the future universe evolu- 
tion to finite-time singularity [HI [7] . The classification of the (four) finite-time future singularities 
has been made in Ref.|B]. Some of these types future singularities are softer than other and not all 
physical quantities (scale factor, effective energy density and pressure) necessarly diverge at this 
finite future time. 

The presence of finite-time future singularities may cause serious problems in the black holes 
or stellar astrophysics [9]. Thus, it is of some interest to understand if any natural scenario to cure 
such singularities exists. 

In the present paper, we generalize the results obtained in Ref. [6_ and [7J, and we explore the 
role of viscous fluids [TT] within this class of modified f{R, G) gravity models, investigating how 
the singularities may change or disappear, due to the contribution of these viscous fluids. 

The paper is organized as follows. In Sec. II, we present the model of f(R, G)-gravity and the 
equation of state for generic viscous fluid. Hence, we show the gravitational field equations and 
the energy conservation law of fluid. In Sec. Ill, we explain the four types of the finite-time future 
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singularities. Next, in Sec. IV, we study the behaviour of fluids with constant thermodinamical 
parameter u on finite-time future singularities. Moreover, in Sec. V, we explore the singularities 
in the presence of a generic class of fluids with u non constant. In Sec. V, we consider quantum 
effects related to singularities in the realistic Hu-Sawicki Model of /(i?)-modified gravity. Finally, 
conclusions are given in Sec. VI, and in the Appendix we discuss the energy conditions for wich 
the singularities appear. 

We use units of &b = c = H = 1 and denote the gravitational constant Gm- 



II Formalism 

We consider the simple action of f(R, G)-gravity: 

~R + f(R,G) 



S= d 



■74, 



16ttG 



C 



N 



(ILl) 



where g is the determinant of the metric tensor g^ u , and C is the fluid Lagrangian, and f(R, G) is 
a function of the Ricci scalar R and the Gauss-Bonnet invariant G = R 2 — AR tlv R ,tv + R llv ^ a R ,lv ^ a 
(R^ u and R^^ are the Ricci tensor and the Riemann tensor, respectively). The general flat FRW 
space-time is described by the metric: 

ds 2 = -dt 2 +a 2 (t)dx 2 , (II.2) 

where a(t) is the scale factor of the universe. The Hubble parameter is H = a/a (the point denotes 
the derivative with respect to the cosmic time t) and for the above metric, the scalar curvature 
and the Gauss-Bonnet read: 

R = 6 (2H 2 + H^j , (II.3) 

2 ( rj2 



G = 24FT [H z + H) . (II.4) 

The scalar expansion is V^u M = 3H , being the four- velocity of the cosmic fluid introduced 
in Ea. pi.ip . We assume the fluid equation of state in the form: 

p = u(p)p - 3HC( P ) , (II.5) 

where p and p are the pressure and energy density of fluid, respectively. The thermodynamical 
variable uj(p) is an arbitrary function of the density p. (,(p) is the bulk viscosity and in general 
it depends on p. On thermodynamical grounds, in order to have the positive sign of the entropy 
change in an irreversible process, has to be a positive quantity, so we assume > 0. For 
the stress-energy tensor T^ Vl one has : 

+ {u{p)p - 3HC(p)) (fiv + Uu,u v ) . (II.6) 
The equations of motion (EOM) are derived from Eq. (|H.ip by using the variational principle [S]: 

PG + P = ^G- N H ^ W 
PG+P= -7T7T- (iH + 3H 2 ) . (ITS) 



8ttG 



The part of modified gravity is formally included into the modified energy density pa and the 
modified pressure pa as follows: 

Pg = -^^{ 24 # 3 /g + SH 2 f R + 6Hf' R + (/ - Rf R - Gf'a)) , (II.9) 
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-i— {8ff 2 / G + 2/^ + 4F/^ + A6Hf' G {H + H 2 ) + /^(4H + 6H 2 ) + 

(f-Rf' R -Gf' G )}. (11.10) 

Here, we have used the following expressions: 

_ df(R,G) _ 0/(fl,GQ 

/fl " ai? ' /g - ■ (IL11) 

Note that the effective energy density of the universe is p e // — Pg + P an d the effective pressure 
is Peff = Pg + P- whereas for the effective equation of state, we obtain: 

- Pe ff - i 2H m 191 

We // = ~ ~~ 1 ~TEH- (11.12) 

Pe// 3ii z 

The fluid energy conservation law is a consequence of the EOM (|II.7|) - (|II.8|) : 

p + 3Hp(l + uj(p)) = 9H 2 ((p), (11.13) 
and, in general, we write (^(p) like a function of H, C(H). 

Ill Future-time singularities 

We are interested in the modified f(R, G)-models that produce some type of finite-time future 
singularities (see, for instance, |6J,[7J), giving by the Hubble parameter: 

H =w^ +H ^ (IIL1) 

where h, Hq and to are positive constants and t < to because it should be for expanding universe. 
/3 is a positive constant or a negative non-integer number, so that, when t is close to to, H or some 
derivative of H and therefore the curvature become singular. Hq could be relevant in the limit 
t — > to only when /3 < 0, and we can assume Hq — when f3 > 0. 

Note that such choice of Hubble parameter corresponds to accelerated universe, because on 
the singular solution of Ea. (|III.l[> the strong energy condition is violated (see the Appendix). 

If /3 = 1, Eq. (|III.l[) with Hq = implies that the scale factor a(t) behaves as: 

<*(*) = ? a ° m • (III.2) 

W (to - t) h y ' 

Ifj8>0andj8^1, Ea. lpTTjl with H = yields: 

h{t -t) 1 -' 3 

a(t) = a e . (III.3) 

In the end, if j3 < 0, we have to take into account the constant Hq > in Eq. pil.l[) . and the 
general form of a(t) is given by: 



-( to -*)(ffo- ht y_r g ) 



o(t) =a e V P / . (III.4) 

Here, ao(^ 0) and in what follows po are positive constants. 

The finite-time future singularities can be classified in the following way[8j: 

• Type I (Big Rip): for t — > to, a(t) — > oo, p c ff — >■ °o and \p e s\ — > oo. The case in which p e g 
and p c ff are finite at to is also included. It corresponds to /? = 1 and /3 > 1. 
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Type II (sudden): for t — > to, a(t) — » do, p c ff — > Po an d |p c ff| — > oo. It corresponds to 
-l</3<0. 

Type III: for i — > to, a(t) — > ao, /O c ff — >• oo and |p c ff | - ^ oo. It corresponds to < < 1. 

Type IV: for f — >• £o> a (0 — ► ao, Pcff —> 0, |p G ff | — > and higher derivatives of H diverge. The 
case in which p and/or p tend to finite values is also included. It corresponds to < —1 but 
is not any integer number. 

We note that in the present paper, we call singularities for = 1 and those for > 1 as the 
"Big Rip" singularities and the "Type I" singularities, respectively. 

In the absence of fluids, the effective density and pressure of the universe are given by modified 
gravity, so that p e ff — pa and p e ff = Pg- A f(R, G)-model shows some types of singularity, if 
Ea. (|III.lj l is a solution of the EOM, namely: 

3 f h 2 2hH al 

PG = ^XJ^W + K^W + i ' ( } 

1 f 2fe/3 3fe 2 6feg OIj2l 

PG = ^8^1 (^Tjm + + + 3H } ■ (IH.6) 

We want to see how the presence of viscous fluids influences the behaviour of singular f(R, G)- 
models (i.e. models that in absence of fluids produce some singularities and satisfy Ea. (|III.5jl - 
(IIII.6jO . We will check the solutions of the fluid energy density when H is singular by using the 
energy conservation law of Ea. (|II.13|) . and we will see how changes the effective density (and, as 
a consequence, the effective pressure) of the universe, and if the singularities are still realized. In 
particular, we are interested in the quintessence (—1 < u) < 0) and phantom (uj < —1) region 
(referring to an ideal fluid), and we assume that the contribute of ordinary matter and radiation 
in expanding universe is too small with respect to the modified gravity. 

We investigate the cases of uj constant and uj dependent on energy density p. 



IV uj constant 

Let us start considering the simple case when w is a constant. We take different choices of bulk 
viscosity £. 

A. Non-viscous case 

In the non- viscous case £ = (perfect fluid), the solution of Eq. (|II.13l) assumes the classical 
form: 

p = p a- 3 ^ , (IV.l) 

where po is a positive constant and a is the scale factor. As a consequence, on the sigular form of 
H of Eq.lUnj}, P behaves as (see Ea. (jlTL^ - (jlTL^ - (|lTL4l) ): 

P = Po(to - tf h( - 1+ ^ , when = 1, (IV.2) 

P = Poe ^ , when j3 > , /3 ^ 1 , (IV.3) 



o 3(i+u<)(fo-t) (go- M Vi ~ ) 



p-poe V J 5 when^<0. (IV.4) 

For = 1 (Big Rip) and > 1 (Type I singularity), p grows up and becomes relevant when 
t is close to to only if u) < — 1. It means that phantom fluids increase the effective density and 
pressure of the universe in the case of Big Rip and Type I singularities, whereas quintessence fluids 
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(uj > — 1) become negligible and do not influence the asymptotic behaviour of f(R, G) models that 
realize this kind of singularities. 

In Einstein's gravity (f(R,G) = 0), Ea. pVT2j) and Ea. pT7j) admit the solution: 



H = - - 7 -, (IV.5) 

3(l + w)(t -t) 

and we can see that the phantom fluid produces the Big Rip for H = h' /(to — t), where h' = 
-2/3(1 + w). 

In general, in f(R, G)-modified gravity, in the presence of phantom fluid, the asymptotically 
Big Rip singularity could appear if pq in Eq. pi.9p diverges less than H 2 (~ (to — t) 2 ) on the 
singular solution of Eq. pV.5l) , namely the modified gravity becomes negligible with respect to the 
fluid energy density in Eq. pi.7p . On the other hand, if a f(R, G) model realizes the Big Rip for 
a certain value of h, the fluid energy density p of Eo. pV.2p becomes negligible on this singular 
solution if ui > —(1 + 2/3h), because it diverges less than H 2 in Ea. pi.7p . 

When f3 > 1, the energy density p of phantom fluid diverges exponentially in Ea. (IIV.3p . so that 
the EOM (|II.7[) - flll.8|) become inconsistent and the Type I singularity is not realized in f(R,G)- 
gravity. 

When < j3 < 1, p tends to po with time in Ea. (IIV.3p . and it is asymptotically negligible with 
respect to H 2 (~ (to — t)~ 2/3 ). In this case, a f(R, G)-model realizing Type III singularity, is not 
influenced by perfect fluids on this kind of singularity. 

For Type II and IV singular models (/3 < 0) , the presence of quintessence or phantom fluids can 
make the singularities worse. Note that H 2 of Type II and IV singularities tends to the constant 
3Hq /8ttGn like ~ (to — t)~^, while pin Eq. pV.4p . after the developing of the exponential function 
in power series, tends to po like ~ (1 + ^)(to — t). 

In the case of— 1 < /3 < 0, a large value of energy density po becomes relevant in Eo. pi.7p and 
could change the numerical value of Ho for which the singularity appears in f(R, G)-gravity, but 
does not necessarly avoid the singularity. 

In the case of (3 < —1, the dynamical behaviour of Eo. pi.7p could become inconsistent, be- 
cause p behaves as (to — t) and it is larger than the time-dependent part of H 2 (~ (to — In 
particular, the Type IV singularities with \/3\ >> 1 are very difficult to realize in the presence of 
phantom or quintessence fluids. 



Examples: 

• In the model f(G) = — a\fG 1 where a is a positive constant, the Type I singularity or the 
Big Rip for some values of h > 1 could occur [7]. If we add a phantom fluid (u < —1), the 
Type I singularity is avoided, while the Big Rip could still appear. 

If u> < — 5/3 (namely, u < — (1 + 2/3/i) for any value of h > 1), the fluid energy density of 
Eq. pV.2p grows up faster than H 2 when h > 1, and this kind of Big Rip is not realized. On 
the other hand, the phantom fluid could produce the Big Rip for some value of < h' < 1, 
when h! = —2/3(1 + cj) like in Eo. pV.5p . However, it is possible to verify, by using Ea. pi.9p . 
that pc of this model, when < h! < 1, diverges still like H 2 , but is negative, so that, if its 
modulus is larger than the fluid density, the effective energy density of the universe becomes 
negative and the Big Rip is avoided, because Eq. (III.7[) results inconsistent. It depends on 
the value of a parameter with respect to the fluid density po. 



• The model f(R) — uR? , where a is a constant, could realize a Type II singularity when 
7 < or a Type IV singularity when 7 > 2 [B]. In both cases we assume Hq negligible in 
Ea. llITLTl) . 

The presence of quintessence or phantom fluids does not avoid the Type II singularity, 
because the numerical value of Hq changes on the singular solution (Hq = ^8ttGnPo/3), 
but the dynamical behaviour of the modified function f(R) keeps the same, due to the fact 
that R tends to infinity in Ea. pi.3p . and is not influenced by the costant H - Moreover, if 
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we use a phantom fluid, there is the possibility that the Type II singularity is changed into 
the Big Rip in the form of Eo. pV.5|) . because, when H ~ (to — t) -1 , it is easy to verify that 
Pq of the model tends to zero, so that the fluid is dominant and makes the future singularity 
stronger. 

When 7 > 2, the Type IV singularity could be avoided by phantom or quintessence fluids, 
especially if 7 is very close to two. This is because the model is singular for (3 = 7/(2 — 7), so 
that |/3| >> 1 if 7 is close to two. As a consequence, other future scenarios for the universe 
are possible. For example, if 7 = 3, the model admits an instable de Sitter solution with 
Rds — -\/l/ a |16J, or the phantom fluid may produce an accelerating phase. 

• The model /(G) = -aG 7 , where a > and 7 > 1, shows the Type II singularity with Hq 
negligible and —1 < (3 < —1/3 in Eq. (|III.l[) . Now, the presence of phantom or quintessence 
fluids with large density po , avoids the Type II singularity. The value of Hq and the dynamical 
behaviour of /(G) change together, because in the case of Ho = 0, when — 1</3<— 1/3, G 
tends to zero in Eg. ([II. 4ft . but if Ho is substantially different to zero, G diverges to infinitive 
and Eq. (III.7[ ) for this kind of model becomes inconsistent on the Type II singularity. 



B. Constant viscosity 

Now, we introduce bulk viscosity in cosmic fluid. Note that viscous fluids belong to more general 
inhomogeneous EoS fluids introduced in Ref.|10|. 

Suppose to have the bulk viscosity equal to a constant, £ = £o- Ea. (|II.13[) yields: 

p = p a- 3{1+ ^ + 9Coa~ 3(1+w) f dt'a(t') 1+3u a(t') 2 . (IV.6) 



For the Big Rip (fi = 1), p behaves as (see Eq. plL^ ): 

, = Mt0 -^,« ) + _i^L_. (IV . 7) 

The bulk viscosity is not asymptotically relevant on the Big Rip solution, because the Hubble 
parameter H 2 diverges like l/(to — t) 2 , while the viscosity part of p diverges more slowly, like 
l/(f - t), and in the EOM pI7 ]) -(|lLg ]l we can neglect it. 

In the paper |11| is written the general exactly form of H(t) for fluids with constant viscosity 
in absence of modified gravity (f(R, G) = 0): 

{ ) ~ o 1 3/1 1 ; /?n — — gggggziy ■ ( IV - 8 ) 

3 + 3 (1 + ^)V24^GVp l 12jGwCo 

H{t) shows a finite-time future singularity when t tends to to, where to = (127rGjv£) -1 ln(l — 
247tGtvCo/(1 + w ) V247rGjv/9o)- If we expand the exponential functions around to, we obtain: 

Hit) ^ - * ... 1 ^ + + O{to t) , (IV.9) 

3(1 + uj) (t - t) l+oj 

that corresponds to Eq. pil.ll) with f3 — 1 (Big Rip), h = —2/(3 + 3w), where uj < —1, and 
Ho = 87tGatCo/(1 + w). The viscosity Co is not relevant in the asymptotic limit of H (here, Hq is 
negative, but the first positive term of H is much larger), and we recover Eq. pV.5|) . that is valid 
for phantom perfect fluids. 

In order to study the effects of the viscosity on Type I, II, III and IV singular models, it is 
worth considering the asymptotic behaviour of the conservation law in Eq. pi,13[) . We require that 
the left part diverges like the right part on the singular solutions: 



it - ty 1 "V " [to - typ + [to - ty 



p + Ml + ^i- — + H ) ^ + - ^+9H 2 ( Q , (IV.10) 



where we take Hq — if /? > 0. In what follows, we neglect the homogeneous solutions, already 
discussed in the previous chapter. 

The asymptotic solutions of Eq. (|IV.10|) arc: 

p-j- ^ when/3>l, (IV.ll) 



(1 + w)(to - t)P 

(2/3-l)(i -t)2^-i ' 
9hH C , 3J? Co 



whenl>/3>0, (IV. 12) 



when P< 0,H ^0. (IV. 13) 



p " (^-^(to-i)^ 1 1 + f 

In the first case (/? > 1), it is possible to see that p diverges more slowly than H 2 , so that 
viscous fluids do not influence the asymptotically behaviour of Type I singular models in Eq. pi.7[ ). 
due to the constant viscosity. 

Also in the second case (0 < j3 < 1), viscous fluids are asymptotically avoidable in the case of 
Type III singular models, because Ea. (|IV.12p diverges less than H 2 . 

In the end, we consider fluids that tend to a non-negligible energy density when (3 < 0. It 
automatically leads to Hq ^ in Eq. pi.7[) and p behaves as in Eq. pV.13[) . Large bulk viscosity 
Co becomes relevant in the EOM, determining the value of Hq for which the singularity occurs in 
some model of f(R, G)-gravity. Moreover, if u < — 1, the effective energy density (namely, pc + p) 
could be negative and avoid the Type II and IV singularities for expanding universe (Hq > 0). 

Example: 

• The model f(R) — a(e~ XR — 1), where a and A are positive constants, shows the Type II 
singularity in the form of Eq. pil.l[) with Hq — and /3 very close to — 1 + [7]. 

A fluid with u> > — 1 and constant viscosity Co large with respect to a, changes the value of 
Hq (Hq ~ -y/8"7rGArCo/(l + UJ ))i but the Type II singularity can still occur. 
On the other hand, a fluid with u> < — 1 and Co >> ct, makes the singularity unphisical (the 
solution of Eq. pi.7jl - pV.13[) leads to Hq < and the singularity is for contracting universe). 

C. Viscosity proportional to H 

This is the case £ = 3Ht. As £ is assumed to be positive, the constant r has to be positive. 
Ea. llrmj) yields: 

p = p Q a- 3il+uj) + 27ra- 3{1+u) J dt'a(t') 3uJ a(t') 3 . (IV.14) 
For the Big Rip (/3 = 0), p behaves as: 

27h 3 r 

P= (tQ-t) 2 (2 + 3h + 3hco) ' (IV ' 15) 

and diverges like H 2 in Eq. pi.7[) . so that, in general, the Big Rip could still appear in the f(R, G)- 
models that realize this kind of singularity. 

In absence of modified gravity (f(R, G) = 0), Eq.jTVT5| and Eq. pi.7p admit the solution: 

2 1 

H = (72ttG n t - 3(1 + oj)) (t - t) ' (IV ' 16) 
and realize the Big Rip for H = h'/(to — t), where h' = 2/ (72ttGnt — 3(1 + uj)). h' is positive 

(1 +lu) - 2AitG n t < 0. (IV.17) 

It means that, in f(R, G)-gravity, it could appear the Big Rip driven by phantom fluid or fluid 
in the quintessence region with sufficiently large bulk viscosity, if the modified gravity on the Big 
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Rip solution is small in Ea. (|II.7[) . namely pc diverges less than H 2 given by Eo. pV.16p . On the 
other hand, if (1+u) — 24-kGnt > 0, the fluid does not realize the Big Rip for expanding universe. 
The other asymptotic solutions of Eq. (|II.13[) are: 



9// 



2„ 



27rh 3 

(3/3-l)(to-t)W-i ' 

TJhElr 9H2t 



when (3>1, (IV. 18) 

when < f3 < 1 , (IV. 19) 



^ (/^jfaV^ TT^ when^O.Fo^O. (IV.20) 

For j3 > 1, p diverges like H 2 if oj > — 1. Thus the fluid could asymptotically produce the Type 
I singularity when the modified gravity is negligible or behaves as p on this kind of solution. On 
the other hand, if uj < —1, for large values of viscosity r, the Type I singularity could be avoided 
in f{R, G)-gravity, because the effective energy density of the universe becomes negative, making 
inconsistent Eq. (|II.7[) . 

When < fj < 1, the fluid does not influence the Type III singular models and can be neglected 
on the singularity, because Ea. (|IV.19p diverges less than H 2 . 

When f3 < and Hg ^ 0, the bulk viscosity can influence the singular solutions for large values 
of r/(l + uj) in Ea. pV.20p . In particular, if cj < —1, the effective energy density of the universe 
could become negative, avoiding Type II and IV singularities in f{R, G)-gravity. 



Examples: 

• The model f(R, G) — —a(G/R), where a is a positive constant, shows the Type I singularity [7]. 
A fluid with uj > —1 and energy density in the form of Eq. (|IV.18p . influences the feature of 
the singularity by changing some numerical value, but the Type I singularity is still realized. 
In addition, if t is sufficiently large, Eq. pV.17p is satisfied and an other possible scenario is 
the Big Rip. 

If uj < — 1, large values of r make negative the effective energy density of the universe on 
the Type I singularity, that could be changed into the Big Rip. 

• In the model f(R,G) = a(G/R), where a is a positive constant, the Type III, II and 
IV singularities could appear[7|. The presence of fluids with ui < — 1 and large viscosity 
proportional to H, does not influence the Type III singularity, but could change the Types 
II and IV into the Big Rip, like in the previous example. 



V uj not a constant 

In this general case, uj = u>(p), we are interested in some simple solution of Eq. (|II.13[> . when 
H is singular and behaves as in Eq. pil.ip . We consider viscous fluid, whose thermodinamical 
parameter uj is given by: 

oj = A p a - 1 -l 1 (V.l) 

where Aq(^ 0) and a are constants. When a — 1, we find the case when w is a constant. Let us 
suppose the following form of bulk viscosity £: 

( = (3H) n r. (V.2) 

Here, n is a natural number and r is a positive constant different to zero. 
The energy conservation law leads: 

p + 3HA Q p a =9H 2 {3H) n r, (V.3) 

from which we may get the (asymptotic) solutions of the fluid energy density when H is singular. 
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In what follows, we consider several examples. 



For the Big Rip singularity (/3 = 1), some simple (asymptotic) solutions of Eq. (|V.3[) are given 

by: 

P = ~, ; „, , , . : — —r for a = 1 , (V.4) 



, on+lun+l 

P - -T7 \~~TT for a > 1 . (V.5) 



Ea. (|V.4p corresponds to the cases when u is a constant. For n = 0, 1, we find Eq. pV.7[) and 
Ea. (|IV.15[) . When a = 1 and n > 1, the fluid energy density diverges faster than H 2 (~ (to — 1)~ 2 ) 
and the EOM (|II.7|) - (|II.8|) become inconsistent on the Big Rip. We can say that fluids with ui 
constant and bulk viscosity proportional to H n , where n > 1, avoid the Big Rip in f(R, G)-gravity. 
The same happens in the presence of viscous fluids with uj non constant and n + 1 > 2a, where 
a > 1, as in Eq. (|V.5p . 

For Type I singularities (j3 > 1), an asymptotic, simple solution of Eo. (|V.3|) is: 



3 n+1 h n+1 T 
4)(*o-t) (n+1)/5 



for a ^ 1 . (V.6) 



The cases a — I and n = 0, 1 correspond to Ea. pV.ll[) and Eq. pV.18p . The fluid avoid the 
Type I singularities if 2a < n + 1 when a ^ 1, so that its energy density diverges faster than H 2 
in Eo. pi.7|) . It means that, if the viscosity behaves as a power function of H larger than one, the 
fluid is able to avoid the Big Rip and Type I singularities in f(R, G)-gravity. 

Note that the viscosity is introduced in the EOM by the fluid pressure of Eq. pi.5p . On the Big 
Rip and Type I singularities, the curvature R behaves as H 2 in Eq. pi.3p . Motivated by fact that 
the correction term R 1 with 7 > 2 cures Big Rip and Type I singularities in f(R, G) gravity [7], 
we may directly conclude that fluid viscosity proportional to H 1 ^ 1 shows the same effect, like we 
have just seen. 

For Type III singularities (0 < f3 < 1), an asymptotic solution of Eq. p/.3p is: 

on+2i,n+2 

forl/2<a<l. (V.7) 



(2/3 + n/3- l)(t - t) 2 ^ 3 - 1 



The cases a = 1 and n = 0,1 correspond to Eo. pV.12p and Ea. pV.19p . The fluid energy 
density diverges faster than H 2 when n > 1//3. In principle, if a f(R, G)-theory shows the Type 
III singularity for a certain value of /3, the presence of a fluid with viscosity proportional to H n , 
where n > 1/(3 (and, as a consequence, n > 1), can make inconsistent the EOM and avoid this 
kind of singularity. Otherwise, it could appear a new Type III singularity realized by fluid for 
H = h/(to — t) 1 / 71 , so that p ~ H 2 , solving in some cases Eq. pi.7p . 

For Type II and IV singularities (/3 < 0), if Hq =^ 0, an asymptotic solution of Eq ^V.3p is given 

by: 



3 n+2 H£ +1 hT f3 n+1 HS +1 T 



(p-l)(to-t)f>->- V A Q 



for a ^ 1 . (V.8) 



The cases a — 1 and n = 0, 1 correspond to Eq. PV.13P and Eq. pV.20p . In general, this kind 
of fluid influences the feature of Type II and IV singularities in /(i?, G)-gravity, but not necessarly 
avoid they. 



If the viscosity is equal to zero (r — 0), Eq. ljV.ip and Eq. (|V.3p yield: 

i 

(a-1) (3A \n^ 
a 



(V.9) 
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where a(t) is, as usual, the scale factor, ao a positive parameter and a/1 (non perfect fluids). 
We may take a > 1 and Aq positive, and consider the parameter ao very small with respect to 
a(t), so that, in general, p is positive. 
In addition, we set 

(3A (a-l))T^=-^_, (V.10) 

07TLTJV 

where Ho is a positive parameter. As a consequence, we have: 

a (ty 



In- 

a 



r 8nG N 

In Einstein's gravity (f(R, G) = 0), the first equation of motion (|II.7|) reads: 

-i 2(a-l)/(2a-l) 



(v.ii) 



a(i) = «oExp|6 



2 

2a-l 



(2a-l)(V3ff t) 
a- 1 



(V.12) 



Note that for large values of a, the fluid energy density tends to H 2 /8ttGn, and Ea. ljV.ip . by 
using Ea. dV.10p . leads to u) ~ —1, and a(£) ~ aoe^ */ 3 ((anti)de-Sitter universe). 

By matching Eq. (|III.3[) and Eq. (|V.12[) , we can see that the fluid produces the Type I singularity 
for j3 = l/(2a — 1), where 1/2 < a < 1, and in general does not avoid Big Rip and Type I 
singularities in f(R, G)-gravity. 

Since on Type III singularity the fluid energy density tends to a constant in Ea. (|V.ll[) . it does 
not avoid this singularity in f(R, G)-gravity, when the effective energy density diverges. 

At last, on Types II and IV singularities, it is easy to see that, by substituting Ea. (|III.4j) into 
the scale factor a(t), the asymptotic time-dependent part of Eq. (|V.ll[) behaves as ~ (to — t) and 
it is larger than the time-dependent part of H 2 (~ (to — t)^ 13 ), when j3 < — 1, in Eq . (|II.7|) . so that, 
in general, especially Type IV singularities become difficoult to realize in f(R, G)-gravity in the 
presence of this kind of fluids. 

A special case of non-viscous fluid with u) non-constant is the Chaplygin gas[13j, whose equation 
of state is: 

P=-^, (V.13) 
P 

where Aq is a positive constant. Ea. (|II.13[) leads: 



P=^o + ^. (V.14) 

Since a(t) does not tend to zero for any type of singularities, the energy density always tends 
to a constant in the limit t — > to, and the Chaplygin gas does not influence the asymptotic be- 
haviour of f(R, G)-models in wich appear Big Rip or Type I and III singularities. In the case of 
f(R, G)-models that realize Type II and IV singularities, the energy density is not negligible and 



could influence the behaviur of the EOM (1II.7I) - (|II.8|> . avoiding Type IV singularities. 
Example: 

• In the model f(R) = aR 1 / 2 , with a positive constant, the Type III singularity for j3 = 1/3 
could appear [7]. A fluid with ui constant and viscosity proportional to H n , where n > 3, 
avoid this kind of singularity. It is interesting to see that in this case, since if /3 = l/n, pa 
diverges faster than H 2 , the fluid does not produce a new Type III singularity, due to the 
contribute of modified gravity. Moreover, the model is free of any other type of singularity 
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VI Quantum effects 

Consider next the quantum contribution to the conformal anomaly. The complete fluid energy 
density is: 

Ptot=P + PA, (VI. 1) 

where pa is given by quantum effects. Taking the trace of the conformal anomaly energy- 
momentum tensor, 

T A = -PA + 3 PA , (VI.2) 
plus observing the energy conservation law, 

PA + 3H(p A +pa) = 0, (VI.3) 

we find that: 

PA = ~ PA ~ m ■ (VL4) 

Thus we obtain for the conformal anomaly energy density [17 : 

PA = -^4 J dta 4 HT a = 

-i J dta 4 H^-12bH 2 + 2Ab'{-H 2 + H 2 H + H 4 ) - 

(46 + 66") (ff + 7HH + AH 2 + 12H 2 H)} . (VI.5) 

Here, b, b' and b" are constants, occuring in the expression for the conformal trace anomaly: 

T A = b(F + ^DR) + b'G + b" R . (VI.6) 

Here, F is the squared Weyl Tensor and G the Gauss-Bonnet invariant. Explicitly, if there are 
N scalars, N±/ 2 spinors, N 2 gravitons and Njjd higher derivative conformal scalars, one has for b 
and b' the following expressions: 

. N + 6N 1/2 + 12N 1 +611N 2 -8N HD 

b = 120(4^ ' (VL?) 

N + llN 1/2 + 62N X + U11N 2 - 28N HD 

b m^T 2 ' ( } 

whereas b" is an arbitrary constant whose value depends on the regularization. 
The quantum corrected FRW Eg. (III. 71) is: 

PG+P + PA^^^H 2 . (VI.9) 

07TG-JV 

Quantum effects become relevant for large values of curvature R and when the effective energy 
density of the universe is not too much large. In particular, this is the case of Type II singularities. 

Let us consider a simple example of realistic /(i?)-model. The Hu-Sawicki Model[14J reproduces 
the current acceleration of universe: 

f(m _ m 2 Cl (R/m 2 ) n _ to 2 ci m 2 ci/c 2 fVTinl 

J[ ' c 2 (i?/m 2 )» + l c 2 c 2 (R/m 2 ) n + 1 ' 1 ' ' 

where m 2 is a mass scale, c\ and c 2 are positive parameters and n a natural number. 
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This is a "one-step model" (for reviews, see [151 116| ) and it is characterized by the existence of 
one transition scalar curvature, so that, for R = the modified gravity vanishes and we find the 
Minkowski solution of Special Relativity, and, for large value of R, the model exibits a constant 
asymptotic behaviour, namely f(R >> to 2 ) ~ — cim 2 /c2- The model has been very carefully 
constructed and c\rv? jc-i ~ 2A, where A is the Cosmological Constant. 

The Hu-Sawiki Model could become singular. In particular, it shows a Type II singularity 
when H behaves as: 

h 

{t~ty 

By solving the asymptotic limit of equations of motion (|III.5[> - (|III.6|) in the absence of fluids, 
when j3 = —n/n + 2, we find: 



H = u _ + H a , -l</3<0. (VI.11) 



6n 2 (n + l) (2 + 



{n + 2) 2 \-Qn 





n+2 

(VI.12) 



(VI.13) 



Here, h is positive if n is an even number. In this case the model could show the Type II 
singularity. Note that Hq = H e is, where Hds is the constant Hubble parameter in the de Sitter 
universe, because the contribute to the constant derives from the term in f(R) into Eq. (|II.9[) . and 
f(R) tends to the Cosmological Constant for large positive or negative values of R. Not necessarly 
a one step model tends to the Cosmological Constant for negative values of R (see, for instance, 
the exponential models in Ref.[TS]). 

In order to avoid the singualrity, we can introduce, for example, a quintessence fluid in the 
cosmological scenario. On the Type II singularity, the quantum effects of fluid become relevant 
and Ea. (|VI.5|) behaves as: 

pa ^w^- (VL14) 

Here, a is a number. This term is very large with respect to the density of quintessence fluid, and 
diverges in Ea. (|VI.9[) . so that the model is protected by singularities. 

Note that, in the de Sitter universe, quantum effects are negligible and the de Sitter solution is 
asymptotically stable, due to the fact that the quintessence fluid disappears in expanding universe. 



VII Conclusion 

In the present paper, we have investigated the finite-time future singularities in f(G, i?)-gravity 
in the presence of viscous fluids. We have considered different fluids and the possibility to avoid 
the singularities. 

In principle, perfect fluids may add singularities in f(R, G)-models or make some singularities 
stronger, but, due to introdution of bulk viscosity, all the four types of singularity could be removed 
from theory. Note that viscous fluid does not mean the modification of gravity. This is just some 
effective fluid added to gravity under consideration and it suggest a possible scenario to protect 
the theory from singularities. 

It is interesting to compare these results with Dark Matter effects to future singularities as 
described in Ref.[18], where it was found that coupling of singular forms of Dark Energy with 
Dark Matter may prevent only from Type II and Type IV future singularities. The viscosity leads 
to high curvature correction terms as a function of Hubble parameter and the effect of viscous 
fluid maybe more similar to R 2 term which may cure all types of singularities as it is first shown 
in Ref.[19]. In this respect, it would be of interest to study the combined effect of Dark Matter 
and viscous fluid to future singularities in modified gravity. 

In addition, we have discussed a possible way to resolve the singularity problem, by taking into 
account quantum gravity effects, that become important in the case of Type II singularities. 
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Appendix 

A Energy conditions near the singularities 

We briefly discuss the energy conditions maybe related with occurrence of singularities. We 
have four types of energy conditions: 

• Weak energy condition (WEC): p > and p + p ^ 0; 

• Strong energy condition (SEC): p + p ^ and p + 3p ^ 0; 

• Null energy condition (NEC): p + p ^ 0; 

• Dominant energy condition (DEC): p ^ \p\. 

On the singular solution H = h/(to — t) 13 + Ho, we have 

I hp . , 

P + P= -**G N (to-*)'* 1 ' (L1) 

where p and p are the effective energy density and pressure of the universe (deriving from modified 
gravity and fluid). 

The effective DE related with Type I and III singularities violate the NEC, whereas DE related 
with Types II and III satisfy the NEC. 
Note that 

3 / , „ hH h 2 hp \ 

The effective DE related with Type II singularities (— 1 < f3 < 0) violate the SEC for small 
value of t, but, when t is close to to, the last term of Eq. (|I.2[) is dominant and the SEC is satisfied 
on the singular solution. 

Finally, in the case of Type IV singularities, when t is really close to to, the term H 2 could be 
dominant and the SEC is violated, expecially if |/3| << 1. 

At last, it is easy to see that, on the singular solutions, when t is near to to, the DEC is always 
violated except for large value of Ho in the case of Type IV singularities, but also in this case the 
behaviour of universe approaching the singular solution violate the DEC. 

We note that in general viscous fluid curing the singularities not necessarly conflict with these 
energy conditions on singular solutions. 
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